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Abstract 



For over a century, x-ray scattering has been the most powerful tool for determining the 
equilibrium structure of crystalline materials. Deviations from perfect periodicity, for ex- 
ample due to thermal motion of the atoms, reduces the intensity of the Bragg peaks as well 
as produces structure in the diffuse scattering background [1^ 2j. Analysis of the thermal 
diffuse scattering (TDS) had been used to determine interatomic force constants and phonon 
dispersion in relatively simple cases [3-6j before inelastic neutron scattering became the pre- 
ferred technique to study lattice dynamics [7j . With the advent of intense synchrotron x-ray 
sources, there was a renewed interest in TDS for measuring phonon dispersion [8j. The 
relatively short x-ray pulses emanating from these sources also enables the measurement 
of phonon dynamics in the time domain. Prior experiments on nonequilibrium phonons 
were either limited by time-resolution [9l [10] and/or to relatively long wavelength excita- 
tions [IT-TS]. Here we present the first images of nonequilibrium phonons throughout the 
Brillouin zone in photoexcited III-V semiconductors, indium-phosphide (InP) and indium- 
antimonide (InSb), using picosecond time-resolved diffuse scattering. In each case, we find 
that the lattice remain out of equilibrium for several hundred picoseconds up to nanoseconds 
after laser excitation. The non-equilibrium population is dominated by transverse acoustic 
phonons which in InP are directed along high-symmetry directions. The results have wide 
implications for the detailed study of electron-phonon and phonon-phonon coupling in solids. 

Experiments were performed at the BioCARS beamline of the Advanced Photon Source 
at Argonne National Laboratory. A dual undulator setup provided an x-ray photon flux 
at the sample of ^ 10^° photons/pulse with photon energy of 12 — 15 keV in a < 4 % 
bandwidth. Single 100 ps x-ray pulses were isolated by a high-speed chopper at a rate of 
40 Hz. A Ti: sapphire regenerative amplifier synchronized to the accelerator RF pumped an 
optical parametric amplifier to produce visible pulses 1.2 ps in duration. The pulse energy at 
the sample was ^ 62 /xJ focused on a 0.14 x 2.1 mm spot. The samples were single crystals 
of InP and InSb with the (001) direction perpendicular to the wafer surface. The crystals 
were oriented at grazing incidence so that the x-rays were incident primarily along the (100) 
direction. Typical incidence angles oi a = 0.3^ were used to match the laser excitation depth 
with the volume probed by the x-rays. A large area detector (MarCCD 165 mm diameter) 
placed 64 mm behind the sample collects the x-rays with scattering angles of up to 2^ ^ 51°. 
In choosing the x-ray energy and crystal orientation, special care was taken to avoid Bragg 
reflections on the detector. 
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Figure[T](b) shows a static TDS pattern from InP at room temperature taken with 15 keV 
x-ray photons. The upper half of the image is shown, the lower half being blocked by the 
sample. The large area detector captures a solid angle of the Ewald sphere in reciprocal 
space projected onto a plane. Each pixel has associated a different scattering vector Q 
that depends on the x-ray energy and the sample-detector distance. The intensity pattern 
originates primarily from (inelastic) scattering from thermally activated phonons and shows 
structure that reflects the details of the phonon dispersion in the material and the phonon 
population at a given temperature [14]. 

Ultrafast laser excitation of the material initially produces a non-equilibrium state where 
phonons may not be populated according to the Bose factor. This transient state eventually 
leads to a new equilibrium where the lattice is left at a higher temperature. These pertur- 
bations of the phonon population modify the equilibrium diflFuse scattering pattern shown 
in Fig. 1 (b). As a function of laser-x-ray delay, t^, we record a frame with laser exposure, 
/(tj), and a dark frame, /(off). Figures [l] (rf) — (g) show the change in the image upon laser 
excitation for several representative time delays [i.e. I{tj) — /(off)]. Each data frame I{tj) 
is obtained by averaging 100 individual images with 100 single-pulse exposures per image. 
Initially, the areas near the zone-centers brighten at t = ps (d) while the scattering keeps 
increasing several hundred ps after excitation (e) and (f), before decreasing after a few ns 
(g), eventually returning to equilibrium. The laser and x-rays are coincident at t = ps (d), 
however the initial increase in scattering takes ^ 100 ps to develop (e). 

We consider the intensity of one-phonon thermal diffuse scattering given byp^ 

where ujj{q) is the frequency of the phonon mode in branch j with reduced wavevector q, 
nj(q) is the phonon population and Fj{Q) is given by 

FjiQ) = E -7=^~'''(^ ■ (2) 

s y'^s 

In this expression, /g, m^, Ms are the atomic scattering factor, the mass, and the Debye- 
Waller factor [15j of atom s at position r^, e^j^q is the phonon polarization vector, and 
Kg is the closest reciprocal lattice vector to Q, i.e. Q = q + Kg. In Eq. ([T]) we omitted 
additional slowly varying factors from the polarization dependence of the scattering as well 
as geometrical factors [14j . 

3 



In thermal equilibrium, nj{q) = -c^^.(q))fcgT -^ ^ which, for the low frequency limit gives 
nj(q) oc ksT /ujj{c\). This multiplied by the l/ujj(c\) factor in Eq. ([T]) gives a strong 
contribution from acoustic modes near the center of the Brillouin zone where their frequency 
goes to zero. These low-frequency phonons give rise to the four brightest diffuse spots in 
Fig. [l] (b). Directions in reciprocal space where the dispersion relation is soft appear in 
Fig. [l] (b) as bright lines connecting the near-zone-center spots. The quantity i^j(Q) in Eq. 
([2]) has the form of a structure factor modified by the dot product (Q • e^j^q), which selects 
phonon modes polarized along the scattering vector. A slight deviation from the four-fold 
symmetry of the crystal is expected between the vertical and horizontal directions due to 
the horizontally polarized x-rays. The eight small spots of a few pixels in diameter in Fig. [T] 
(b) are crystal surface truncation rods which appear due to the grazing incidence geometry 
and are well known from surface diflFraction [16j. 

Equations ([T]) and ^ can be used to compute the TDS if the phonon frequencies and 
eigenvectors are known. The inverse problem, i.e. obtaining and e^j^q from /(Q) 

can be done model- free only in simple cases [17] . To gain insight on the different phonon 
contributions we implemented a Born- von Karman model of the lattice dynamics with force 
constants up to six nearest neighbors [lU [19]. InP and InSb crystallize in the zincblende 
structure with two atoms per unit cell. Six phonon branches and eigen-displacements, three 
acoustic and three optical, are obtained. This model ignores the ionic component of the 
forces and thus does not reproduce the phonon dispersion with all detail, particularly it 
does not give a splitting of the LO and TO phonons at zone-center. However the insight 
gained is extremely valuable for understanding the static diffuse images particularly in sep- 
arating the contributions from transverse (TA) and longitudinal (LA) acoustic phonons. 
Figure [l](c) shows a calculated TDS intensity pattern using this force constant model. This 
calculation agrees well with the static measurements and reproduces the main features. More 
importantly, as we show later, it proves useful in understanding the time-resolved data. 

To separate the different time-scales we performed singular value decomposition (SVD) of 
the time-resolved data. This method has been applied successfully to a wide set of problems 
from small angle scattering [20j and the analysis of spectroscopic data [21 J to time-resolved 
Laue diffraction of protein crystals [22j. The procedure is as follows. We write each difference 
frame I{tj) — /(off) as a column vector Xj and construct the matrix X whose columns are 
the vectors Xj. Thus, adjacent columns of X correspond to consecutive frames in the time- 
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delay series. The SVD states that X = USV^, where S = diag(5i) is the diagonal matrix 
of singular values with 5^ > 0, and U and V are orthogonal matrices whose columns, Uj 
and Vj, are the left and right singular vectors, respectively. The vectors Uj represent time- 
independent "populations", and form a complete orthonormal basis. The vectors vj contain 
the time-dependent information of each of the populations and also form an orthonormal 
basis. The singular values {si} act as weight factors for these vectors and represent the 
amount of "signal" contained in each singular vector in the reconstructed data. 

In Fig. [2] we summarize the results of the SVD analysis for both InP [(a) and (b)] and 
InSb [(c) and (d)]. The truncation rods were masked out to avoid artifacts from saturated 
pixels. Due to the larger lattice constant, the data for InSb were taken at 13 keV to match 
the same reciprocal space covererage on the detector as InP at 15 keV. On the left we 
show the (time-independent) left singular vectors (ISV) for the two most significant singular 
values (SV) with the corresponding (time-dependent) right singular vectors (rSV) scaled by 
the corresponding SV on the right. For both materials, much of the signal is contained in 
the first SV [Fig. [2] (a) and (c)], the time dependence of which consists of a sharp increase 
in the diffuse scattering that decays in few ns time scale (Fig. [2] (a) and (c) right panels). 
It is worth noting here the clear differences between InP and InSb, particularly the lack of 
scattering intensity in InSb near the top center of the first ISV in Fig. [2] (c), which is also 
present in the dark frame (not shown). This is due to the small structure factor for the 
(204) refiection in InSb, which is not small in InP due to the large difference in the atomic 
scattering factors. The shape of the image from this ISV resembles the equilibrium image, 
which suggests this ISV corresponds to lattice heating and diffusion that decays within a 
time-scale of few ns. 

The second singular vectors however show a more complex behavior both in reciprocal 
space as well as in the time domain. In the case of InP, the ISV [Fig. [2](b) (left)] shows sharp 
bright and dark areas indicating that parts of the image brighten while other areas decrease 
their scattering. The fact that the rSV [Fig. [2] (b) (right)] peaks at t = 400 ps means 
that the image on the left continues increasing after the initial excitation which indicates 
that non-equilibrium transient dynamics lead to a delayed heating and cooling of different 
parts of the Brillouin-zone. The same delayed behavior is observed in InSb as seen in Fig. 
[2] (d) (right), which shows a peak at ^ 4 ns. We attribute this behavior to relaxation of 
non-equilibrium phonon populations. However, note that from the way the population nj{q) 
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and frequency (jJj{q) enter in Eq. ([T]), it is difficult to distinguish a population increase from 
a softening of the phonon frequencies. 

Although the (equilibrium) TDS distributions for InP and InSb are similar, as well as 
the first singular vector in Fig. [2] (a) and (c), which refiects similarities in the phonon 
dispersion, the difference in the second ISV is quite remarkable and cannot be explained 
by differences in the equilibrium phonon dispersion alone. It is instructive to separate the 
contribution from different phonon branches to the diffuse scattering. Each phonon branch 
contributes one term to Eq. ([T]), which we compute using the Born- von Karman model 
described above. The top and bottom rows of Fig. [3] show a zoomed view of the second 
ISV in Fig. [2](b) and (d), respectively. Superimposed, we show contour plots that represent 
the calculated contribution to Eq. ([T]) from the LA [(a) and (b)] and the TA [(c) and (d)] 
branches separately. As stated earlier, the largest contribution to the scattering comes from 
areas that are close to the center of the Brillouin-zone, where the contours lines are denser. 
This is further confirmed by superimposing the images with the corresponding Brillouin-zone 
boundaries as shown in (e) and (f). 

From Fig. |3] we see that in InP the areas with strong contribution from TA phonons 
[Fig. |3] (c)] appear bright in the ISV, while areas with more contribution from the LA 
branch [Fig. [3] (a)] are dark. This ISV has the time-dependence shown in the right panel of 
Fig. [2] (b), which means that even at time-delays t > 400 ps TA phonons are still populating 
while the LA branch decreases. In InSb the situation is different. The LA branch shows no 
appreciable decrease in scattering above the noise level [Fig. [3] (b)], while the TA branches 
populate delayed [Fig. [3] (d)]. To further identify which phonon modes are involved within 
the Brillouin-zone, we sketch in Fig. [3] (g) the Miller indices of selected Brillouin-zones for 
InSb at 13 keV. Note that at this energy the number and size of Brillouin zones covered by 
the detector in InSb is similar to InP at 15 keV. Using the diagram in Fig. [3] (g) we can 
assign the InP TA phonons in Fig. |3] (c) to modes along high-symmetry (111) and (010) 
or equivalent directions. InSb on the other hand shows a less directional scattering from 
the TA branches, indicating that the population is distributed more isotropically within the 
Brillouin-zone. We attribute this difference to the large gap between the optical and acoustic 
modes in InP, which leads to less available phonon decay channels [23j. 

In conclusion, we have shown that time-resolved non-thermal x-ray diffuse scattering 
yields heretofore inaccessible information of the non-equilibrium dynamics of solids under 
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ultrafast excitation. InP and InSb show complex non-equilibrium redistributions of the 
scattered intensity throughout the Brillouin-zone that are clearly separated by singular value 
decomposition. This non-equilibrium lattice persists up to several hundred ps. This behavior 
is interpreted using a Born-von Karman model which reveals that TA phonons populate 
several hundred ps after initial laser excitation and, in the case of InP, is very directional 
and appears concomitant with a decrease in the scattering by LA phonons. In the current 
study, we are limited to dynamics that occur on more than 100 ps time-scale due to the 
x-ray pulse duration. We note that measurements of the initial emission of hot phonons by 
the photoexcited carriers and their subsequent anharmonic decay will become possible with 
new x-ray sources such as the Linac Coherent Light Source at SLAG. As these new sources 
of ultra-bright femtosecond x-ray pulses come online, x-ray diffuse scattering becomes a 
promissing approach to explore the ultrafast dynamics of solids under laser excitation at the 
femtosecond time-scale. 
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synchronized with the x-ray pulse, (b) room temperature TDS image of InP at 15 keV oriented 
with (100) parallel to the x-ray direction, (c) Calculated TDS from a Born-von-Karman model with 
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FIG. 2: Results of the singular value decomposition analysis of the time-resolved data. The 
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FIG. 3: Contribution to the TDS from TA and LA branches calculated using Eq. Q. The images 
correspond to an expanded view of the second eigenvectors in Fig. [2] (b) (top row) and Fig. [2] (d) 
(bottom row). The contours in (a) and (b) show the calculated contribution from LA phonons in 
InP and InSb, respectively, (c) and (d) show the contribution from TA phonons for InP and InSb, 
respectively, (e) and (f) shows the same images superimposed with the Brillouin-zone boundaries 
for InP and InSb. (g) Miller indices for selected Brillouin-zones of InSb. 
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